arXiv:l504.04312vl [cond-mat.soft] 16 Apr 2015 


Molecular simulation study of polar order in orthogonal bent core smectic liquid 

crystals 

Stavros D. Peroukidis, Alexandros G. Vanakaras, and Demetri J. Photinos 
Department of Materials Science, University of Patras, Patras, 26504, Greece 

We explore the phase behavior and structure of orthogonal smectic liquid crystals consisting of 
bent-core molecules (BCMs) by means of Monte Carlo molecular simulations. A simple athermal 
molecular model is introduced that describes the basic features of the BCMs. Phase transitions 
between uniaxial and biaxial (antiferroelectric) orthogonal smectics are obtained. The results indi¬ 
cate the presence of local in-plane polar correlations in the uniaxial smectic phase. The macroscopic 
uniaxial-biaxial transformation is rationalized in terms of local polar correlations giving rise to polar 
domains. The size of these polar domains grows larger under the action of an external vector field 
and their internal ordering is enhanced, leading to field-induced biaxial order-disorder transitions. 


I. INTRODUCTION 

Bent-core molecules (BCMs) are known [IHl] to self 
organise into a variety of soft complex structures, often 
showing hierarchal ordering on different length scales. Of 
particular interest for electro-optical [5] and non-linear 
optical devices [5] is the polar (ferroelectric or antifer¬ 
roelectric) switching of nematic and smectic BCM liquid 
crystals (LCs). The tilted stacking of the BCMs in layers 
is rather common and is the defining characteristic of the 
tilted smectic (SmC) class of materials as well as of the 
so called tilted cybotactic nematic class (Ncybc)) com¬ 
bining the presence local tilted molecular layering with 
macroscopic positional disorder [7]. The tilted ordering, 
however, complicates the device architecture, due to the 
monoclinic (or lower) symmetry it confers to the materi¬ 
als. In addition, it gives rise to the formation of chevron 
structures and defects which can be a serious disadvan¬ 
tage for the use of tilted LC materials in electro-optic 
devices. On the other hand, when the stacking of the 
molecules into layers is not tilted, dehning the so called 
orthogonal LCs, whether of the smectic (SmA) class or 
the orthogonal cybotactic nematic (Ncy^), the afore¬ 
mentioned tilt-inflicted drawbacks are removed and such 
materials could be very advantageous if they show polar 
order and optical biaxiality [5]. 

The conventional nomenclature for orthogonal smectic 
phases of BCM materials includes the SmA, the SmAPi^, 
the SmAP^, the SmAP^jy^ and the SmAPo,. The SmAPy 4 
was first observed by Eremin et. al. [9] and characterized 
as an orthogonal biaxial smectic phase exhibiting anti¬ 
ferroelectric switching behaviour and consisting of polar 
layers of antiferroelectric order. The SmAP^j as well as 
the SmAPfla were studied by Pociecha et. al. [10] who 
proposed that the polar directors of the layers are ar¬ 
ranged randomly in the first case and are helicoidally 
modulated along the layer normal, in the second. A sim¬ 
ilar interpretation for the structure of SmAP^j has been 
provided by Panarin et. al. m- A different model of the 
structrure of SmAP^j has been proposed by Shimbo et. 
al. |12| wherein the molecules are organized in polar clus¬ 
ters (domains) that are randomly distributed within the 
layers, giving rise to a macroscopically uniaxial and ap- 


olar state. Furthermore, the recently observed SmAP^A 
is suggested to consist of randomly aligned nonpolar do¬ 
mains of local antiferroelectric ordering [T3| . 

These different interpretations of the connection be¬ 
tween local structure and the macroscopic behavior of 
BCM LCs have stimulated a lot of interest for theoretical 
research. Systematic attempts have been made to ratio¬ 
nalize this behaviour using theory and computer simula¬ 
tions UMI], but there is still much to be done, especially 
on the local structure-macroscopic properties connection. 
In previous works, the BC molecules were modelled as an 
assembly of interacting sites in a bent configuration by 
joining i) two hard spherocylinders |16] . ii) Lenard-Jones 
spheres iii) Gay-Berne particles m, and iv) soft 
spherocyliders m- A fully atomistic simulation of BCM 
nematics has also been reported [20]. 

In the present work, we have introduced a hard-core 
model for the molecular interactions that permits us to 
model BCMs with relatively short arms and sharp bend 
angles. The phase behavior and local structure is ratio¬ 
nalized in terms of the shape anisotropy together with 
athermal specific interactions [22] • Finally, the response 
to external stimuli is also discussed and compared with 
experimental findings. 


II. MOLECULAR MODEL AND SIMULATION 
DETAILS 

We have employed a coarse grain molecular model 
which captures two essential features of the BCMs; these 
are: (i) the anisotropic shape of the mesogenic bent core 
and (ii) the chemical differentiation between the aromatic 
bent core of the mesogen and its terminal alkyl chains. 
The bent core of the molecules (see Fig. 1) consists of two 
spherocylinders (c-segments), of aspect ratio L* = L/D, 
jointed rigidly at their end caps to form the bend an¬ 
gle 7 , where D is the diameter of the cylinders and the 
hemispherical caps. Two terminal spherical segments (t- 
segments) of radius R* = R/D are attached tangentially 
at the free ends of the c-segments of the molecule. In this 
work we have focused on systems with L* = 2, R* = 0.7 
and 7 = 120° in accordance with the size and shape of 
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common bent core smectics. We have also investigated 
systems with different L*, R* and 7 to examine the in¬ 
fluence of these parameters on the phase behaviour. The 
chemical incompatibility of different molecular units has 
been modelled by introducing the following differentia¬ 
tions among the segmental interactions [ 22 ]: (i) the c- 
segments interact with all types of segments (c or t) via 
hard core repulsion on overlap and null interaction other¬ 
wise, and (ii) the interaction potential of a t-t pair of seg¬ 
ments is null, irrespectively of mutual overlapping. Our 
results demonstrate directly that these minimal molecu¬ 
lar features are sufficient to generate the experimentally 
observed polymorphism and phase sequences. 

The phase behaviour and molecular organization of the 
so defined BCM systems is studied by means of Metropo¬ 
lis Monte Carlo simulations in the isobaric isothermal 
ensemble (NpT) using variable size simulation boxes 
with periodic boundary conditions |23|. We have per¬ 
formed compression and expansion series by systemati¬ 
cally varying the pressure. Equilibration requires on the 
order of 2x10® cycles and a further 5xl0®-lxl0® cycles 
to be used for the calculation of ensemble averages of 
quantities of interest. A MC cycle consists, on aver¬ 
age, of N trial attempts (translations, orientations and 
translations-orientations) of a randomly chosen particle 
and one volume change attempt. The temperature and 
pressure are expressed in reduced units. The reduced 
temperature is kept constant T* = ksT/e = 1, where T 
is the temperature, fcs is the Boltzmann constant and 
e is the unit of energy. The reduced pressure is given 
by P* = pD^ /ksT. The number density is defined as 
p* = ND^/V, where N is the number of molecules and 
V is the volume. The inter-molecular distances r and cor¬ 
relation lengths ^ are scaled by the diameter D, defining 
the reduced distances r* = r/D, the correlation lengths 
= ^/D and the respective reduced scattering vector 
magnitudes q* = 27r/r*. 


III. RESULTS AND DISCUSSION 

A. Phase behaviour and molecular organisation in 
the smectic phases 


Systems consisting of BCMs with L* = 2, R* = 0.7 and 
7 = 120° (see Fig. 1), comprising from N = 845 to 3364 
molecules have been simulated. To quantify the orienta¬ 
tional order and to identify the principal axes frame of the 
simulated systems we have diagonalized the order-tensors 


where 


[IS|b Q“ = 2N Sill 3 (“» ■ -®) “ 

a = {x,y,z) represents any of the molecular axes and 
A,B = X,Y,Z denote the axes of the simulation 
box. The eigenvector associated with the largest posi¬ 
tive eigenvalue of the three order-tensors is taken as the 
primary director n. 

At relatively low pressures, the order parameters 
S'" = ^2 i&i ■ n)) approximately vanish, which 


indicates the absence of orientational order. Here 
P2 (di • n) denotes the 2nd Legendre polynomial of the 
direction of the x or y or z molecular axis relative to 
the primary director. The biaxial order parameters are 
evaluated from the components of the tensor = 



di 7 ^ bi are axes of the ith molecule and I, m denote the 
secondary principal axes of the phase. 
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FIG. 1. Coarse-grain molecular model of BCM used in the 
simulations. The molecular axes frame is shown. 



In Fig. 2(a) we present the calculated pressure vs num¬ 
ber density p* — p* equation of state (EoS) and in Fig. 
2 (b) the the pressure dependence of the order parameters 
for a system oi N = 845 molecules. The system under¬ 
goes a transition, accompanied by a density jump, from 
the isotropic (I) to an orthogonal smectic (SmA) phase. 
A characteristic snapshot of this smectic phase is shown 
in Fig. 3a, from which it is evident that the molecular 
axis 2 is ordered. The respective order parameter in¬ 
creases from nearly zero to approximately 0.9 (see Fig. 
2(b)). At even higher pressures a transition to a biax¬ 
ial smectic (SmAH^) with antiferroelectric arrangement 
occurs. A characteristic snapshot is shown in Fig. 3(b). 
The biaxial order parameter of the SmAF^ is in the range 
0.4 < < 0.8 (see Fig. 2(b)). The reverse process, 

i.e. expansion from the SmAH^i that has been obtained 
by compression, shows considerable hysteresis, with ex¬ 
tended range of stability for the orthogonal smectic SmA 
phase, compared to the compression series. Finally, at 
lower pressures the SmA is destabilized in favour of the 
isotropic phase. 


The structure of the isotropic liquid and of the ordered 
phases is examined through: (i) the calculated two di¬ 
mensional x-ray scattering pattern |24j ; the intermolecu- 
lar scattering intensity ImteriQ) reported in this study, 
is given by subtracting the single molecular intensity 
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FIG. 2. (a) The equation of state for systems consisting of 

N =845 BCMs. (b) The order parameters and as a 
function of the reduced pressure p* . Squares and circles cor¬ 
respond to expansion series from ordered phases and triangles 
to compression series from the isotropic phase. 



FIG. 3. Representative snapshots for the simulated BCM 
systems: (a) SmA at p* = 1.1 and (b) SmTPA at p* = 1.6. 


Imoi{q) from the total intensity Itot{q), Iinter{q)= Itot{q) 


- imoiiq) with itot{q) = fT{q)fTiq) and imoliq) = 
Fj (q), where /t(q) is the total scattering (struc¬ 
ture) factor ^3 (q) ' '*'1 )> (q) is th® spher¬ 

ical uniform scattering factor of the jth molecule, Vj is 
its position and (ii) a set of positional two-dimensional 
pair correlation densities m defines as: 

ffo n) ^ ^ (ra - Vij ■ a,) S (ra - r,j ■ b,^ 


V 


V 2 I 

[r^3 ■ 

1 Oi xbij 




which give the molecular density on the plane defined 
by the axes of a single molecule. Here 0(x) denotes the 
step-function (0 = 1 for x > 0 and 0 = 0 otherwise), 
and cr = Zl/2. 

First we study the structure of the thermodynamically 
stable isotropic phase (away from the hysteresis regime) 
which lacks long range positional and orientational or¬ 
der. The ( 7 g’^ {y, z) and (x, z) show broad maxima 
(depicted by red arrows in Fig. 4) which are located at 
a distance slightly lower than one molecular length. The 
observation of these maxima indicates positional correla¬ 
tions along the z molecular axis. A side by side packing 
of the molecular cores is evident from the intense maxima 
indicated by green arrows in Fig. 4(a). 

It is interesting to note that, within the range of 
stability of the isotropic phase, the location of the 
maxima on the pattern does not change on increasing 
the pressure; only their intensity increases. The polar 
intermolecular correlations are rather weak and rapidly 
decay to zero with distance; this is evident from Fig. 
4(c), showing the calculated orientational correlation 


function {x,y), defined according to g'^'^ {ra,rb) = 

( (** -^j)^ (»'a - rij ■ a^) 6 (rb- ■ 5,^ 
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Nevertheless, the nearest neighbours of a molecule in 
its X — y plane appear to point in the same direction. 
The molecular arrangement described by the correlation 
functions in real space is also supported by the calculated 
x-ray scattering pattern. A spherical scatterer is assigned 
to the molecular apex (i.e. at the origin of the molecular 
axis frame). The inner ring (small angle region) corre¬ 
sponds to interlayer distance r* = 27r/q* w 6.0 in accor¬ 
dance with the intermolecular distance calculated from 
the density functions. Finally, the outer ring corresponds 
to the side-by-side molecular packing. In the isotropic 
phase the correlations are local, excluding any long range 
orientational or positional order. This, to our knowl¬ 
edge, is the first local structure analysis in the isotropic 
phase of BCMs systems using molecular simulations. In¬ 
terestingly, the same smectic-like clusters, i.e domains of 
enhanced positional correlations, have been observed in 
real systems of bent core molecules in the isotropic phase 
[25] and were invoked for the interpretation of interest¬ 
ing properties such as large flow birefringence [50]. At 
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higher pressures the system undergoes a transition to an 
orthogonal smectic phase. This phase consists of apolar 
layers with spacing approximately equal to that observed 
in the local positional correlations. 





FIG. 4. Calculated correlation functions for systems of 
BCMs with L* = 2, i?* = 0.7 and 7 = 120° in the isotropic 
phase at p* — 0.7: (a) {x, z), (b) [y, z) in real space, 

(c) {x, y) and (d) x-ray scattering pattern with scattering 

vector magnitude q* = 211 /r*. 


In order to examine the in-plane polar order of the 
smectic phase we have evaluated the gi'^ {x, y) functions 
(see Fig. 5), from which we have found that short range 
polar correlations are present and that the range of ori¬ 
entational correlations is strongly anisotropic (extending 
in the x direction more than in the y). Accordingly, each 
molecule in the layer can be viewed as being surrounded 
in the x — y plane by an anisometric region contain¬ 
ing neighbouring molecules that have biaxial correlations 
with that molecule. We have estimated the polar corre¬ 
lation lengths along the x and y molecular axes 

from the g^'^ {x,y) functions as the optimal parameters 
for fitting the separate x and y dependences of this func¬ 
tion to the functional forms g^'^ {x,0) = exp{—\x\/^x) 
and gi'^ (0, y) = exp {—\y\/^y). The correlation lengths 
Q and increase with increasing pressure. These are 
roughly equal with very small values (^ l.OI?) just above 
the isotropic-smectic transition and increase to = 4.2 
and ^y = 2.5 at pressures deep in the SmA phase (see Fig. 
6 ). The calculated correlation lengths suggest that each 
polar domain incorporates up to a few tens of molecules. 

Representative snapshots of the directions of the arrow 
vectors (identified with the unit vectors of the molecular 
X axis in Fig. 1) within a smectic layer are shown in Figs. 


7(a)-(b). There it can be seen that, in the SmA phase, 
small polar domains exist and are distributed randomly 
within the layer. Some polar domains are indicated in 
Fig. 7(a)-(b) using ellipsoids with major and minor ra¬ 
dius of approximately the lengths and , respectively. 
The fluidity of particles and the isotropic positional dis¬ 
tribution within the layers should be noted. Optical in¬ 
spection of uncorrelated snapshots (every one thousand 
of MC cycles) show a continuous formation and disap¬ 
pearance of polar domains. This indicates that these 
domains can not be considered as stable floating ’’polar 
islands” within the layers. The stochastic nature of the 
Monte Carlo simulations do not allow inferences about 
the dynamics of the domain formation and disappear¬ 
ance process. The size of these domains increases with 
pressure (this is also reflected on the persistence of the 
{x,y) function over longer intermolecular distances). 
Interestingly, the sign alternation of (x, z) with z 
near the transition to the SmAF^ indicates that the small 
polar domains positioned one just above the other on 
successive layers show antiferroelectric order correlations 
(see Fig. 5(b)). Therefore, the SmA phase in this limit 
incorporates the characteristics of a conventional SmA 
phase and also consists of weakly correlated antiferro- 
electrically ordered domains. We term this smectic state 
as SmApc. Note that, according to ref. 113: the lay¬ 
ers in a SmAFp are polar, with the polar directors ran¬ 
domly arranged across the different smectic layers, thus 
producing a macroscopically apolar smectic phase. This 
interpretation differs from the one given in ref. nana, 
which is supported by our findings, i.e. each of the layers 
is overall apolar and consists of randomly oriented polar 
domains. Hence, the origin of the macroscopic apolarity 
of the phase is attributed to cancellations of local polar 
correlations (extending few molecular diameters) within 
each layer, as opposed to the mechanism of polarity can¬ 
cellations among macroscopically polar layers. 



FIG. 5. Gontour maps of calculated correlation functions for 
systems of BGMs with L* = 2, A* = 0.7 and 7 = 120°, (a) 
{x,y), (b) g^'^ {x,z) for the SmApc phase at p* = 1.30. 


Further compression of the system brings it to a transi¬ 
tion from the SmApc to the SmAFyi phase. In this phase. 
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FIG. 6 . Calculated polar correlation lengths Q (squares) and 
Q (circles) vs pressure in the SmApc phase. Solid symbols: 
without external field, open symbol: under external field a* = 
0.09. 
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each layer becomes macroscopically polar, with the direc¬ 
tion of polarity alternating from layer to layer, thus form¬ 
ing a macroscopic antiferroelectric structure. This im¬ 
plies that the polar domains observed in Sm4 Pc increase 
in size and merge upon the transition to the Sm4P^ giv¬ 
ing rise to long-range two-dimensional polar order within 
each layer. In this case the antiferroelectric coupling be¬ 
tween adjacent smectic layers counterbalances the spon¬ 
taneous tendency towards disorder in two-dimensional 
systems. We note here that the phase sequence: I- 
SmT-Sm^P/j-Sm^P^-Cr has recently been verihed ex¬ 
perimentally in our case, the Smj4-Sm4Pfj phase 
transition of [27] appears as a nearly continuous transi¬ 
tion from an orthogonal smectic with small polar domains 
to an orthogonal smectic with larger ones. 


On increasing the bend angle to 7 = 140° or 160° for 
the same L* = 2 and R* = 0.7, the range of stability of 
the Sm4 phase increases. On the other hand, increas¬ 
ing the length of the arms to L* = 3 with R* = 0.7 
and 7 = 120° destabilizes the Sm4 phase, which is now 
observed within a small range of pressures. Finally, a 
slight decrease of the size of the t-segments destabilizes 
the liquid crystalline phases: the system with L* = 2 
and R* = 0.5 and 7 = 120° exhibits only an isotropic to 
crystal transition. 


B. Domain structure in the presence of static 
alignment polar field 

A macroscopic measure of the size and internal order¬ 
ing of the clusters can be provided by the susceptibil¬ 
ity of the system to field-induced ordering |5S1 A 
similar measure is provided by the threshold value of 



FIG. 7. Representative snapshots presenting the molecular 
arrow vector (unit vector along the molecular x axis) in a 
single smectic layer at two pressures for systems with (c,d) 
and without (a,b) external aligning field, top panel, a* = 
0: (a) SmApc phase at p* = 1.20 and (b) at p* = 1.30; 
bottom panel: same pressures with a* — 0.09. The ellipses 
are drown to indicate some polar domains of the systems and 
to emphasise the larger polar correlation lengths in the field- 
on systems. 


the applied field that is necessary to produce a tran¬ 
sition from the state of randomly oriented clusters to 
the macroscopically ordered phase m- The transition 
mechanism in this case entails the merging and size in¬ 
crease of the clusters. For clusters identified on the basis 
of a vector property, such as the polar ordering of the 
molecules considered in these simulations, the appropri¬ 
ate susceptibility refers to the direct coupling of a vec¬ 
tor property reflecting molecular polarity with a vector 
field. For electro-optic applications the response of the 
LC phase to the electric field is of primary interest and 
therefore the relevant coupling would be that of a per¬ 
manent molecular dipole moment with an applied electric 
field. However, as the aim of the present simulations is 
to identify the possible phase organization modes that 
originate directly from the bent shape of the molecules, 
the model molecules are not endowed with permanent 
electrostatic moments or polarisabilities. Strictly, there¬ 
fore, the only dipole moments of the molecules are the 
steric ones, associated with their bent shape. In principle 
these would couple with mechanical vector fields produc¬ 
ing, for example, polar flow-alignment. Alternatively, it 
may be formally assumed that the molecules carry per¬ 
manent electric dipole moments that couple directly to 
an externally applied electric field but are very weak to 
generate signihcant intermolecular interactions. Gener¬ 
ally, electric dipole interactions are known to influence 
the phase behavior and the local structure of LC sys¬ 
tems [3T]. Recently, a SmA-SmAPi? phase transition has 
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been predicted theoretically [H] in a system of bent-core 
molecules possessing a permanent dipole moment along 
the arrow direction. An extension of the present molec¬ 
ular simulations is under way for the study of the phase 
behavior of systems of bent-core molecules possessing a 
permanent dipole, and of their response to an external 
electric field. 

In the present work we investigate the response of 
the system to an externally applied polar-field E that 
couples linearly to the molecular arrow vector, x, and 
contributes to the potential energy of the ith molecule 

as Ui = —a (^i ■ , where a is a coupling parame¬ 

ter. The dimensionless strength of the field is mea¬ 
sured as a* = a/kT. When the field is off, the sys¬ 
tem does not show a net spontaneous polarity. This 
is clearly reflected through the vanishing first rank po¬ 
lar order parameter Pf = Here m 

denotes the secondary principal axis of the phase along 
which the molecular arrow-vectors order in the presence 
of the external field (a* 0); it coincides with the di¬ 

rection of the applied field. The dependence of biaxiality 
and polar order parameters on pressure, for various field 
strengths, are presented in Figs. 8(a) and 8(b) respec¬ 
tively. From these plots is easily concluded that fields 
with strengths a* > 0.06 render the low-pressure SmApc 
macroscopically polar and therefore biaxial. In contrast, 
a measurable macroscopic polarity in the SmAP^i phase 
is obtained only when the applied field exceeds a critical 
strength a* « 0.12. Fields with strengths up to a* = 0.12 
do not induce any detectable polarity in the crystalline 
phase {p* > 1.9). We stress here that the field strengths 
with a* < 0.12 leave the calculated macroscopic biax¬ 
ial order parameter of the SmAP^ phase practically un¬ 
changed with respect to the field-off calculations, see Fig. 
8 (a). In addition these fields do not affect the average 
density of the systems. These observations suggest that 
the chosen strengths, 0.06 < a* < 0.12, correspond to 
relative weak fields since they do not influence the main 
order parameters of the unperturbed (field-off) systems. 


Application of fields with strengths in the range 0.06 < 
a* < 0.12 to the SmApc phase have the following conse¬ 
quences: i) the size of the polar domains, as quantified 
through the polar correlation lengths and increase 
with the strength of the field, see Figs. 6 and 7(c)-(d), ii) 
the development of orientational correlations between the 
polar domains become long-ranged, leading to a net po¬ 
larity and biaxiality and iii) the polar susceptibility of the 
system increases substantially as the system approaches 
the SmApc-SmAPyi phase transition at p* = 1.30. 

Up to this pressure -where the SmApc-SmAPyi phase 
transition occurs for the purely hardcore system with¬ 
out external field- the macroscopic polarity of the filed- 
on system, as reflected on Pf, grows continuously and 
reaches its maximum value. Above this pressure, depend¬ 
ing on the strength of the applied field, we observe either 



* 

P 


FIG. 8. (a) Biaxiality vs pressure and (b) Polar order pa¬ 

rameter vs pressure for various values of the field strength a* 
for a system of A = 3364 molecules. 


an abrupt drop of the polarity (for a* < 0.09) or, for 
higher fields, a continuous growth of Pf with a tendency 
to saturate to its maximum allowed value as the systems 
reaches its crystallisation pressure. Clearly a* « 0.12 is 
the lowest required field strength to have a field induced 
transition form an orthogonal anti-ferroectically organ¬ 
ised SmA phase (SmAP^) into a ferroelectric orthogonal 
smectic (SmAP^). Fields bellow this threshold leave the 
SmAP^ phase practically unperturbed as indicated by 
the vanishing magnitude of the polar order parameter. 
At even higher pressures {p* > 1.90) an antiferroelec- 
tric crystal is observed. In this case, much higher field 
strengths are expected for the formation of ferroelectric 
crystal state. 

Interestingly, the field response obtained in these sim¬ 
ulation for the SmAP phase is in qualitative agreement 
to the experimentally observed m response to an ex¬ 
ternal electric field for the SmA-SmAPj^ region. Our 
simulations suggest that the underlying microscopic pic¬ 
ture of this transition involves polarly correlated in-plane 
domains that change their size and align under the influ¬ 
ence of the external field, as opposed to randomly ori¬ 
ented macroscopically polar layers which are aligned by 
the external field. 
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IV. CONCLUSIONS 


The thermodynamic stability and the local structure 
of liquid crystalline phases formed by model bent core 
molecules have been examined by the MC-NpT simula¬ 
tion technique. The minimal representation of the inter- 
molecular potential, through selective short-range ather- 
mal interactions, allows us to examine molecular struc¬ 
tures with relative sharp angles and short arms with¬ 
out using the more complex parameterization of previous 
studies [IIHID] • The molecular symmetry of the rigid bent 
core model we have used facilitates the investigation of 
orthogonal smectic phases. The simple molecular model 
reproduces qualitatively experimental observations on or¬ 


thogonal bent core systems m- We have observed the 
presence of local polar in plane domains which are ran¬ 
domly distributed in the SmTpc and their size increases 
on approaching the transition to the Sm^P^. Accord¬ 
ingly, the SmApc-SmAP^ transition immerges from the 
growth of the size of these domains at the transition. 
The response of the uniaxial smectic phase to an exter¬ 
nal vector field that couples to the polar ordering of the 
molecules is sensitively influenced by the presence of po¬ 
lar in-plane domains; these increase their size and get 
aligned by the field. Our results elucidate conflicting in¬ 
terpretations [I1I17] regarding the microscopic mecha¬ 
nisms that are responsible for the macroscopic behavior 
of orthogonal bent core smectics especially for the SmA- 
SmAPp-SmAP^ transition. 
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